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Bimetric gravity theories allow for many different types of cosmological solutions, but not all of 
them are theoretically allowed. In this work we discuss the conditions to satisfy the Higuchi bound 
and to avoid gradient instabilities in the scalar sector at the linear level. We find that in expanding 
universes the ratio of the scale factors of the reference and observable metric has to increase at 
all times. This automatically implies a ghost-free helicity-2 and helicity-0 sector and enforces a 
phantom dark energy. Furthermore, the condition for the absence of gradient instabilities in the 
scalar sector will be analyzed. Finally, we discuss whether cosmological solutions can exist, including 
exotic evolutions like bouncing cosmologies, in which both the Higuchi ghost and scalar instabilities 
are absent at all times. 


I. INTRODUCTION 


The question whether the graviton can have a mass has been asked for a long time and its answer has always been 
accompanied by uncertainties. The linear theory of a massive gravity was first analyzed by Fierz and Pauli (Tj. Since 
then the van Dam-Veltman-Zakharov discontinuity [2] [3] and the appearance of the Boulware-Deser (BD) ghost [1] 
has been challenging the theory. Recently, a theory of a massive spin-2 field was presented in which the coupling 
between an additional fixed tensor field and the metric has a specific structure and is free of the BD ghost EH32 
(see Refs. jl3i 114| for recent reviews on massive gravity). To promote this theory of a massive gravity to a bimetric 
theory, Hassan and Rosen considered a dynamical tensor field where its kinetic term has the same Einstein-Hilbert 
structure as g and does not introduce the BD ghost [ITJ [T5j . This bimetric theory is described by the action 


S = - 


4 

* [ d A xy/=gR{g) - * f d A xy/^fR(f) + f d^xy/^g^ 0 n e n (X) + f d^xyfXjC 

J J n—0 J 


(i) 


where we already set the Planck mass for / M „ to M g (see Refs. |lfil fTT| for further explanation), absorbed m, the 
mass scale of the graviton, into /3 n and expressed masses in units of M g . The interaction between both tensor fields is 
determined by the elementary symmetric polynomials e n of the eigenvalues of the matrices X“ = yjg a P fpj, multiplied 
by arbitrary real coupling constants 0 n . It is convenient to express these free parameters in units of the present Hubble 
expansion rate, Hq. 

A remarkable property of bimetric gravity theories is the possibility of nonstandard, self-accelerating cosmological 
solutions and the ability of making predictions that are different from ACDM. Some of these might be useful for future 
measurements in order to distinguish standard ACDM from bigravity. To benefit from that, one has to pay the price 
and needs to disentangle all the nonviable models from the viable ones. 

Even though this theory has five free parameters, it is not clear whether viable models exist (except for /3i = ... = 
04 = 0 which is simply ACDM) and, if they do, what they look like. In Ref. [T§] simple criteria of viability were 
considered and viable background solutions were presented (see also Ref. HU)- One choice of the coupling parameters, 
in the following simple model, will usually lead to several different cosmological solutions [16l il8li22] . In the following, 
every possible solution will be called a branch. We distinguish between different types of branches, depending on how 
the ratio of the scale factors r of the metrics and gy, v evolves. In solutions on finite branches the ratio evolves 
from zero towards a finite asymptotic value, whereas on infinite branches r becomes infinitely large at early times 
and decreases with time. We call all other branches exotic branches, these usually describe bouncing cosmologies or 
a static universe in the asymptotic past or future. 

So far, only finite and infinite branches were studied in the literature. While many of these are in good agreement 
with observational data at the background level naunms], most of them suffer from scalar instabilities |2lH?7] . It 
seems that only one specific class of models, the infinite-branch bigravity (IBB), is free of scalar instabilities |24| . 
These models are specific infinite branch solutions in which 02 and 03 vanish. Moreover, IBB agrees very well with 


Electronic address: koennig@thphys.uni-heidelberg.de 




2 


observations at the background and linear level [S3J HU. Unfortunately, the authors in Ref. [29, noted that the Higuchi 
bound is generally violated in the early time limit. This bound, first derived in Ref. m, ensures a healthy helicity-0 
mode of the graviton. A violation leads to the appearance of the Higuchi ghost, named after Higuchi who found that 
a spin-2 particle with mass m and 0 < m 2 < 2 H 2 in a de Sitter space leads to a negative norm ESEU (see also Ref. 
[23j in which the Higuchi bound was derived for arbitrary spatially flat FLRW metrics in massive gravity). Note that 
even though IBB seems to be well behaved at the linear level, the appearance of the Higuchi ghost may only be visible 
at higher orders or maybe even only in the full solution |34| . Furthermore, it was found that cosmological solutions 
on this infinite branch suffer from a ghost in the helicity-2 sector at early times [35]. 

The analysis of viable backgrounds in Ref. m : that leaded e.g. to the exclusion of solutions on the exotic branch 
or a vanishing @2 and (j$ in the infinite branch, are, however, based on assumptions like the existence of a matter 
dominated past or the absence of poles in r ', where the prime indicates the derivative with respect to the e-folding 
time t. Even though it would be probably difficult to get exotic solutions in agreement with observations, they are 
a priori not excluded. Moreover, poles in r' = could have a very physical meaning: If r' reaches a pole, then dt 
becomes zero and the Universe undergoes a bounce. Such an example model is shown in Fig. [i] 1 . 

In the following, we will first briefly discuss the background evolution in Sec. [TT| before we then analyze conditions 
for the absence of the Higuchi ghost and scalar instabilities (Sec. HIpV I to draw conclusions about the viability of 
all theoretically possible solutions (Sec. [V|. 




Figure 1: Example of a model (/3i = (0, 0.3, —0.8, 1, —1)) that describes a bouncing universe. Here, the asymptotic past of 
this universe is described by a root at r ~ 0.8. It then contracts, i.e. dt < 0, until r reaches the pole and, finally, expands 
towards a root at r ~ 1.9, which describes a de Sitter point. 


II. EQUATIONS OF MOTION AT BACKGROUND LEVEL 


To find the cosmological background evolution, we 
equations of motion, 


vary the action © with respect to both metrics and find the 


~ \g^R+ (- 1 )"/ 3 " 

n—0 

1 _ 1 ^ 
n—0 

where the overbar denotes curvature of and are suitable polynomials (see Ref. [20] for their definitions). 

At the background level, we will use a Friedmann-Lemaitre-Robertson-Walker (FLRW) ansatz for both metrics with 
two different scale factors, a and b , together with two different time parametrizations t and t = Xt. Throughout this 
work, t represents the e-folding time and a prime denotes the derivative to it. With this ansatz for the metrics, 

g^dxYdx” = a 2 (— H~ 2 dt 2 + dx 2 ) , (4) 

f^dx^dx" = b 2 ( y —X 2 'H~ 2 dt 2 + dx 2 ) , (5) 


9ii\ Y tn) v ( ^9 a/3 fp-tj + 9v\ Y (n)v ( \J 9 aP fp 


=T 

^ UV’) 


( 2 ) 


f^ Y (n)u (y/7°^gp 7 ) + U>y£i)n 


= 0 , 


(3) 


1 Note that this specific model is not viable due to a negative Ti 2 and is only shown for motivation purposes. 


















3 


where 7i is the dimensionless conformal Hubble function, we obtain the goo and /oo equations 

3"H 2 = a 2 {p + fo + 3 for + 3 for 2 + for 3 ) , 

3 'H 2 = a rX 9 (fo + 3for + 3for 2 + for 3 ) . 

{r' + r) 


( 6 ) 

(7) 


Here we introduced the ratio of the scale factors r = b/a. As usual, both the Friedmann and acceleration equations 
for are degenerated with the conservation of the energy, 

p' = -3p{l +w tot ), (8) 

where Wtot denotes the equation of state (EOS) parameter of the total energy density, while there is no extra constraint 
from the acceleration equation for / M „ due to the missing coupling to the energy-momentum tensor. The combination 
of this set of equations leads to 


X = 1 + —. ( 9 ) 

r 

Replacing this constraint in the equations of motion yields 

3 H 2 = a 2 {p + fo + 3for + 3for 2 + for 3 ) , (10) 

a 2 

3'H 2 = — {fo+ 3 for + 3 for 2 + for 3 ) . (11) 

The second alternative Friedmann equation is particularly interesting since it directly determines the evolution of the 
scale factor if the evolution of r is known. 

The sign of & is a priori unknown and, therefore, r could be negative. However, odd powers of r are always 
proportional to either fo or fo. All cosmological solutions with negative r due to a negative scale factor for f^ are 
therefore equivalent to those with positive r after the redefinition fo n +\ —t —fo n +i- From now on, we will assume 
r > 0 2 . The combination of both Friedmann equations leads to an equation for the density as a function of r only, 

p = fo r' 1 - fo + 3fo + 3 {fo - fo) r + {fo - 3fo) r 2 - for 3 . (12) 

It will be useful to study r' , which can be written as 

r' = — = —3 (1 + Wtot) (13) 

P,r P,r 

where we used Eq. ^ in the last step. 


III. HIGUCHI GHOSTS 

Bimetric theories are called ghost-free since the specific structure of the potential term in the Lagrangian avoids an 
additional degree of freedom (d.o.f.), which usually would be the BD ghost. This, however, does not imply that all 
d.o.f. of the massless and massive graviton are not ghosts. 


A. Higuchi bound 

Bimetric gravity theories describe a mixture of a massless and massive spin-2 field. The latter carries five dofs, 
including one helicity-0 mode. In pure massive gravity around a de Sitter spacetime, Higuchi derived a bound for 
the graviton mass to ensure positive norm states noma, a negative norm would imply a ghost helicity-0 mode 


2 This assumption might only be unjustified if both positive and negative values of r are reached at some time. In the later discussion 
we will find that this requires finite, nonzero values of r' at r = 0 in order to produce viable branches. It turns out that these specific 
models will not be able to produce viable cosmologies. 
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and is usually dubbed an Higuchi ghost. The condition for its absence in bimetric gravity theories around a FLRW 
background was derived in Ref. m . 3 In our notations, the bound is 


- (/3i + 2 /3 2 r + /3 3 r 2 ) (l + r 2 ) > fir + 3/3 2 r + 3/3 3 r 2 + /3 4 r 3 = 3r ( — 


which is equivalent to 


/3i + 3r 2 (/3i - /3 3 ) + 2r 3 (3/3 2 - /3 4 ) + 3r 4 /3 3 > 0. 


Interestingly, using Eqs. (12)-(13 I leads to the simple bound 

P,r < 0. 

This condition for the absence of the Higuchi ghost was already derived in Ref. [3Jj] (see also Ref. EH)- Since 

p = —3 (1 + w tot ) — 

nr *' 


(14) 


(15) 


(16) 


(17) 


and p > 0 together with 1 + Wtot > 0 (we are usually considering a combination of pressureless and relativistic matter), 
the bound is equivalent to 


r' > 0. 


(18) 


Note that this holds even for negative values of r. Therefore, in an expanding universe the ratio of the scale factors b 
and a has to increase at all times in order to satisfy the Higuchi bound. Since r' is negative on all infinite branches 
(18], this directly shows that these branches suffer from the Higuchi ghost at all times and confirms the findings in 
Ref. [25] that the bound is violated at least at early times on infinite branches, i.e. large r. On the other hand, 
all finite branches that produce viable backgrounds are free from the Higuchi ghost since viability in these branches 
enforces r' > 0 m ■ This especially includes the finite branch in the f3\ model, i.e. only /3i ^ 0, which was already 
shown to be free of the ghost in Ref. m- 

The rhs of the bound (141 has to be non-negative at all times. Since we already concluded that r > 0 is a valid 
assumption without loss of generality, the Higuchi bound enforces 


B 2 = Pi + 2/3 2 r + /3 3 r 2 > 0, 


(19) 


where B 2 is simply the derivative of p mg , the modified part in the Friedmann equation (101, with respect to r. 
Therefore, the Higuchi bound is related to the change of the amount of dark energy in our Universe with time. 


B. Phantom dark energy 


It is often useful to study the equation of state parameter (EOS), w mg , i.e. the ratio between the pressure and the 
density, from contributions of the modification of gravity. If we know how the matter density in our Universe evolves, 
then the knowledge of w mg enables us to draw conclusions about the acceleration and even the future of our Universe. 
In Ref. [1 8| we showed that Eq. (191 is directly related to the EOS via 


w n 


= _i _ B^ r , 
Pmg 


( 20 ) 


If Pmg > 0 (which, as observations indicate, should hold at least around present time), then the Higuchi bound enforces 
a phantom dark energy. Every cosmological solution in bimetric gravity should therefore have either a Higuchi ghost 
or a phantom dark energy. 

The property of being a phantom is usually thought to come along with a future instability, the “big rip” [37]. Note, 
however, that the EOS is highly time dependent and tends to —1 in the asymptotic future if it described by a root 
in r ', e.g. in most of the finite branch models. A sufficiently fast increase of w mg could then avoid this instability 
and guarantee a better behaved future. A phantom in bimetric gravity is, therefore, not as frightening as in ACDM. 
Thus, a model implying a phantom dark energy should not automatically be related with a problematic future, much 
less be rejected. 


3 Note that the authors in Ref. m used an overall factor of | in front of the potential term in the Lagrangian, which can be compensated 
for by a redefinition of the /3-couplings. 
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C. Tensor ghosts 

Interestingly, the only factor in the lapse of / M „ that is not strictly positive is r + r'. Thus, the only way to get a 
negative lapse is a negative r '. Therefore, fulfilling the Higuchi bound implies a non vanishing and especially positive 
lapse at all times. 

It was mentioned in Ref. j3S] that the relative factor between the kinetic tensor modes for and / M „ is the lapse 
function of and, therefore, a negative lapse is responsible for a ghost in the helicity-2 sector. We conclude that 
the absence of the Higuchi ghost automatically implies the absence of a ghost in the helicity-2 sector. 

As shown in Ref. [35] > the lapse of directly enters in the friction and mass term of the / Ml/ -tensor perturbation 
equation leading to negative values at early times, which is responsible for a fast grow of the tensor modes |35t [58]. 
This is already a signal of the existence of a ghost. To get such a fast growth in the tensor evolution in accordance 
with observations is a challenging but not undoable task [38]. The main problem, however, is the existence of the 
ghost itself. 


D. Consequences of the existence of ghosts 

A ghost helicity-0 or helicity-2 will have a dramatic impact on the viability of a theory. It will lead to an unbounded 
Hamiltonian from below and allow the existence of particles with positive and negative energies. As expounded in Ref. 
[34) . the vacuum state will immediately decay into positive and negative energy particles. This behavior is enough to 
rule out the underlying theory. 4 It is, therefore, not a question of how problematic the evolution is of a field described 
by the equation of motion. A ghost might influence its evolution in a (more or less) unacceptable way, e.g. through a 
negative friction. However, it is not the possibly ill-behaved solution of the perturbation equations that renders the 
theory unphysical, but rather the absence of a stable vacuum state and interactions with negative energy particles. 
It is even possible that such a system could seem to be completely well behaved at all orders in perturbation theory, 
but the perturbative solution still not converge to the exact solution. An example where perturbation theory is even 
able to hide the negative energy solutions, which are present in the full theory, is discussed in Ref. [34]. 

Since bimetric gravity is only an effective field theory, one might wonder whether a ghost could be harmless in this 
setup or whether a ghost is necessarily excited. This is, unfortunately, not the case. As explained in Ref. [33], only 
modes with positive energy are able to decouple, but not a ghost state since there is no positive energy necessary to 
excite a ghost (see also Ref. [40]). Even in effective field theories (and even if the mass of the ghost lies above the 
cutoff) one has to avoid ghosts at all costs. 


IV. EIGENFREQUENCIES OF SCALAR PERTURBATIONS 

After reducing the number of possible cosmological solutions with the demand of the absence of ghosts, we will 
analyze the behavior of scalar perturbations at the linear level. Even though there are already quite a number of 
works in which similar properties were studied, all these investigations were based on strong assumptions, mostly a 
restriction in the parameter space, fixing the EOS of the matter fluid, or focusing on a specific type of branch. In 
the majority of cases, this is a consequence of the complexity of the perturbation equations. Since the aim of this 
work is to draw conclusions about the viability of the most general cosmological solutions in bimetric theories with a 
FLRW background, we will now work out conditions for the absence of gradient instabilities without resigning from 
generality regarding the parameter space, type of branch and nature, i.e., EOS, of the fluid. 

The set of scalar perturbation equations at the linear level can be reduced to a system of two second-order differential 
equations for two potentials Ej describing the two propagating scalar degrees of freedom |26| (see also Refs. [TH1 ' 221 

Ha ng Earn]), 

5 " + AijE'j + BijEj = 0 , ( 21 ) 

where A l:l and Bij are matrices which depend on the background quantities r, B and the parameters of the models. 
The complexity of this system depends crucially on the choice of the gauge. A very convenient one was used in 
Ref. [ 23 , leaned on Ref. [42]. In this work, we take advantage of the relatively simple 5 perturbation equations that 


4 Note that there are “good ghosts”, e.g. the Faddeev-Popov ghost, which are not related to physical degrees of freedom and are, therefore, 
harmless. 

5 Where “simple” means that printing these equations would fill only a couple of pages. 
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the authors found in this gauge (see Ref. |29| for the derivation and printed equations) and analyze them by using 
the ansatz Hj oc e ut . For simplicity, we assume that the eigenfrequencies ui do not depend on time. This is a valid 
assumption as long as |u//o; 2 | <C 1 holds and was confirmed for all models studied in Ref. p3]. In the subhorizon 
limit, we obtain a surprisingly simple expression for the eigenfrequencies, 


1 \n 


, f (i- 2 + l)(/3i-/3 3 r 2 )r' r 2 (/3i+4 ( 82i'+3/3 3 r 2 ) 


p(«J + l) 


/3 1 +20 2 r+/3 3 r 2 


3 r 3 


- 1 


<■' (-»7r' ^ + 1) (ft - Ar 2 ) - 


3 r 3 


- 1 


( 22 ) 


(23) 


which agrees with all previous, but much more complicated, results for one- and two-parameter models that were 
studied in (24j . As already mentioned in Ref |24| , if we assume dark matter only, then for models in which j3 2 = fa = 0 
this reduces to 


W /3o/3i/34 - 



(24) 


In order to discuss stability, we only need to analyze the sign of w 2 : A negative value would imply oscillating and, 
therefore, stable potentials Ej. If, however, w 2 is positive, then Ej grows quickly with time and even faster as the scales 
become smaller. Such an instability is not compatible with the structure in our Universe and needs to be avoided in 
a viable model. 

Let us now introduce B 2 = + 2(3 2 r + /3 3 r 2 to obtain 


w 


2 


k 2 

3 rp.-.'H 2 




3 rp r 


(25) 


Interestingly, the condition for stability depends on how dark energy (and the density of the cosmic fluid) changes 
but not explicitly on how large it is. We observed a similar property during the analysis of the Higuchi bound. Note 
that B 2 is related to the change of the energy density, p >r , and the Hubble expansion via 


Bo = -- 


1 + r 2 


i rp, r - 2 I — 


(26) 


Together with 


we finally arrive at 


cc 2 = 


— I =r~ i B 2 
r 


B2 ,r 1 


(27) 


2 r' (r (r 2 + l) R 2 ,rP,r - (f ) 2 (3 (r 2 + l) R 2jr + rp, r ) + 6 (f ) 4 ) 
r 2 P,r (rp, r -6(f) 2 ) 


- 1 


(28) 


As we will see later, this expression for the eigenfrequencies will become very convenient when analyzing the stability 
around poles in r / , which e.g. always appear in exotic branches. 

It might be useful to study an expression for lo 2 which does not explicitly depend on the /? parameters but on r 
and its derivatives, like Eq. (241. Finding such an expression is always possible when using a set of five independent 
equations to eliminate all coupling parameters. One possibility is the set of equations for r', r", r'", 7 i 2 and p (note 
that the result will not depend on r'") which yields 


/ k \ 2 a 2 pr 2 (w + 1) [2 (w + 1 )r" + r' (6 w 2 — 2w' + 9w + 3)] — 2 'K 2 r' [r' ((w + 1) (r' — 3™) + rw') — r(w + l)r"] 
\h) 3r(w + l)r' (a 2 pr(w + 1) + 2'H 2 r') 

( 29 ) 
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Here, and in all the following equations, we dropped the subscript in wtot for simplicity. If we are interested in 
analyzing the eigenfrequencies at specific epochs, e.g. radiation dominated era (RDE) and matter dominated era, we 
can assume w — const and obtain 

2 / k \ 2 a 2 pr 2 (w + 1) [2?-" + 3r'(2u> + 1)] + 2'H 2 r l [r (r" + 3wr') — r /2 ] 

yT-LJ 3rr' {a 2 pr(w + 1) + 2'H 2 r l ) 

This leads to the condition 


r' \a 2 pr(w + 1) + 2"H 2 r / ] [a 2 pr 2 (w + 1) (2 r" + (6 w + 3 )r') + 2'H 2 r' (r (r" + 3 wr') — r' 2 )] < 0 (31) 

in order to get stable scalar perturbations, i.e. to 2 < 0. When using the Higuchi bound, r' > 0, the first bracket term 
is always positive and, thus, the second one has to be negative. This is equivalent to 

„ r' 2T~L 2 r' (r' — 3rw) — 3oi 2 pr 2 (w + l)(2w + 1) 

< 2r a 2 pr(w + 1)+'H 2 r' 


where we also used r' > 0. Note that the denominator is always positive. If the numerator would be negative, then the 
bound would especially imply r" < 0. However, this is not generally the case and, thus, the condition for stable scalar 
modes is not automatically equivalent to r" < 0 in contrast to the case for /?o/3i/ 34 models during matter domination 
(see Eq. pi). 


A. Radiation-dominated era 


Even though we will not aim to exclude models which are theoretically allowed but do very likely not reproduce 
observational data (an example would be a nearly static universe that did not have a radiation-dominated epoch), 
it is worthwhile to analyze the conditions when the Universe is filled with either relativistic particles or pressureless 
matter only. 

When radiation dominates, i.e. w — 1/3, the eigenfrequencies simplify to 


co 2 = 


k 2 | a 2 pr 2 (2 r" + 5 r') + 2'H 2 r' [rr" — r' 2 + rr') 


3H 2 r 


\a 2 pr + 2'H 2 r' 


(33) 


In the early Universe, the Hubble expansion is usually driven by radiation, i.e. 37~L 2 — a 2 p. With this approximation, 
the condition for stability in the scalar sector becomes 


> - 


r' (r' + lOr) 
r' + 4 r 


(34) 


For large absolute values of r', which is the case e.g. near a pole, we simply obtain r" > — r'. 

In a previous work [24], we studied the eigenfrequencies for IBB and confined ourselves to a universe filled with 
dark matter only. According to Eq. (24), we concluded stable scalar modes because r' increases with time but stays 


negative until reaching the final de Sitter point. Since r , is always negative in IBB, the condition (341 is not necessarily 
valid anymore. However, we can still use condition (31 ). Here, the product of the first two terms is always positive 
since 


r' (a 2 pr(w + 1) + 2 H 2 r') 


IBB 


9a 2 /3ir (r 2 +1) 


/ (w + 1) (/3i + /3 4 r 3 - 3(3ir 2 ) 
y /3i — 2/3 4 r 3 + 3/Sir 2 


2 

> 0 . 


(35) 


Therefore, we can analyze the third factor and, assuming w £ (—1,1) for simplicity, find that stable modes are 
guaranteed if 


3/3ir 2 </3i+/3 4 r 3 , (36) 

which is not only satisfied in the RDE, i.e. large r (note that both /3i and /3 4 have to be positive in order to get 
a viable cosmological background), but, in fact, is equivalent to the condition p > 0 on that branch and, therefore, 
trivially satisfied at all times. 












B. Matter-dominated era 


Let us study the regime when matter dominates the Universe, 
described through 


Now the EOS vanishes and the scalar modes are 


2 k 2 a 2 pr 2 (2 r" + 3 r’) + 2H 2 r' (rr" — r' 2 ) 

3'H 2 rr' a 2 pr + 2 HPr' 

For stability, we need to satisfy the condition 

„ 2 H 2 r' 3 — 3 a 2 pr 2 r' 

< 2 a 2 pr 2 4 - 2r'H 2 r' 


(37) 


(38) 


If we assume that ILr —> 0 for late times, which should be true when dark energy starts to dominate, then the 
condition of stability reduces to 


r 


// 



(39) 


V. FINDING VIABLE BRANCHES 


We will now raise the question whether branches exist that satisfy both the Higuchi bound and the condition 
for scalar stability. Here we will only focus on cosmological solutions that are not equivalent to ACDM, which of 
course satisfy both conditions. We therefore assume that at least one of the couplings Pi, ...,P± is nonzero. Together 
with conditions of physicality, a, p, 77 2 > 0, we define these as criteria of viability. Note that we allow for solutions 
that have a very nonstandard past, e.g. no matter- or radiation-dominated epoch, or even contracting backgrounds, 
even though these might be hard to compare with observational data. This extends the more restrictive background 
analysis of m- Therefore, not only the finite branch with small r or the infinite one could be viable but also many 
solutions on exotic branches. Many different types of branches exist: some of them start from a root r' = 0, while 
others may evolve from a pole or even pass a pole at some finite time. In many cases it is not directly clear whether 
such branches solve the equations of motion. In particular, every branch always needs to contain a solution of Eq. 
|TT| at present time, i.e. when 77 = a = 1 . 

We start with focusing on finite branches with a root at r = 0. Let us first concentrate on models with Pi 0, 
which always have a root at r = 0 (see Eq. ([13])). In Ref. [24] we generally found scalar instabilities in these type 
of branches. Even though this is based on the assumption of a universe filled with dark matter only, this conclusion 
does not change when considering arbitrary but reasonable EOS parameters. We take the same line of argument and 
study the simple /3i-models, i.e. models with only nonvanishing Pi, since all other models will reduce to these in the 
limit when r gets close to r = 0. The eigenfrequencies in Pi models are given by 


2 1 + 2w — 6r 2 (w + 2) — 9r 4 / k \ 2 1 + 2 w ( k\ 2 

< = -(«) “(3^+ls) (40) 

and, therefore, indicate unstable modes for small values of r as long w > —1/2. Let us consider the previously 
excluded models with Pi = 0 and find 


r —0 


Po ~ 3p2 

h 


(w + 1 ) 


(41) 


Even though the combination Po = 3 P 2 is able to produce a root at r = 0, it will not lead to viable solutions since in 
this case r' = —3 (1 + w) + O (r 2 ) indicates a violation of the Higuchi bound. From this we conclude that 

1. Finite branches with a root at r = 0 always lead to either unstable modes (if Pi ^ 0) or violate the Higuchi 
bound (if pi = 0) for small r. 

On the other hand, r' could be nonzero but still finite at r = 0. In this case, one of the asymptotic points is either 
a pole or the whole branch evolves between two roots at negative and positive r. In the first case, we can assume 
that at least one of the poles is reached at r > 0, otherwise we are able to analyze viability in the “mirrored” model 
corresponding to /tjn+i —> —p 2 n+i- If the branch does not contain any pole, then p^ r has to vanish at r = 0 (roots at 
r/ 0 always indicate a vanishing density whereas a maximum of the density at r 7 ^ 0 leads to a pole). The position 
of the maximum of p at r = 0 requires /3 3 = 0 and leads to p r oc r which cannot be negative for both regions, r > 0 
and r < 0. We can summarize that 
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2. All finite branches with a nonzero and finite r' at r = 0 have to have a pole either in the asymptotic past or 
future. 

Roots, except for those at r = 0, always correspond to a vanishing density. Due to Eq. (13), we will always find a 


pole between two roots ri and r 2 , if both rq and r 2 are nonzero. Therefore, poles could be interesting starting or final 
points of a branch. Whenever such a pole describes the asymptotic future, then r has to go to zero, otherwise the 
pole would not be a stable asymptotic point. Since r' = 'H~ 1 j^r = TL~ 1 ar diverges, we find that 77 needs to vanish 
at this point. On the other hand, if a pole describes the asymptotic past, then we can use the fact that the density 
starts from a finite value. For nonzero val ues of r, this is clear from Eq. (12). It also holds if r = 0 is a pole, since 
this would require fii = @3 = 0 due to Eq. @1 and, therefore, implies pj" 


I r—0 


= 3/32 — Po- If the density is finite at 
early times, the scale factor a has to have aTfhite but nonzero value. In this case, TL needs to be zero at early times, 
too, otherwise one could go backwards in time and we would not have an asymptotic past. Thus, we conclude 

3. Tl has to become zero on a pole, if it describes an asymptotic point. 

Let us assume a pole at r = 0, which, as we already noted, requires /3i = /3 3 = 0 and leads to Tl 2 \ r _ 0 = fi 2 a 2 . From 
the previous conclusion, we need a vanishing 73 2 at r = 0. Note that a > 0, otherwise this would contradict a finite 
density. Therefore, we need /?2 = 0 and, thus, obtain B 2 = 0 for all r, which means that 

4- A pole at r = 0 violates the Higuchi bound. 

For simplicity, we will from now on assume that if there is a pole at r p , then r p > 0. Additionally, we can exclude 
r = 0 from being an asymptotic point due to the previous conclusions. Furthermore, Eqs. (131 and provide the 
limit r' oc —r when taking r — > 00 as long as the density does not vanish (see Ref [18] more detailed explanations). 
This excludes infinite branches, i.e. branches in which r evolves from or to r —> 00 , from being viable due to the 
violation of the Higuchi bound and we find that 

5. The limits r 0 and r —> 00 are no viable asymptotic points. 


We will now consider a root at r ^ 0 as the asymptotic past. Due to Eq. (131, the density vanishes on a root. To 


fulfill the conservation of energy, those models require a contracting universe at early times. If this universe evolves 
to another root (on which again p = 0), then it has to undergo a bounce at p tT = 0 leading to a pole at which 73 = 0. 
Employing the previous conclusions, we find the general statement 

6 . Every viable branch needs to contain at least one pole on which Tl vanishes. 

This result is particularly interesting as it will allow us to draw conclusions when connecting this with the requirement 
of stable scalar perturbations and the absence of the Higuchi ghost. The necessary condition for a pole is p, r —> 0. 
Then, the eigenfrequencies of scalar perturbations around the pole ([28]) reduce to 


W 2 ^ 


XT 2 B. 


2,r“ 


r’ (1 + r 2 ) rp r - 3 (f)' 

r p,*~ 6 (ity 


r 2 p.i 


- 1 


— I vt I 2 B, 


2,r " 


(! 


- 1 


P, r 


(42) 


(43) 


where we used (-y ) 4 <C (-y) 2 , "C (-y)~ (and, additionally, B 2r 7 ^ 0 which, as we will see later, is justified), 

together with 


r Ar-3 (^) 2 _ B 2 (1 + r 2 ) - (^) 2 r 




B 2 (1 + r 2 ) 


^ 1, 


(44) 


which follows from Eq. (26) and B 2 (l + r 2 ) — (^) r ~ i? 2 (l + f 2 ) (note that the Higuchi bound (14) implies 
V 2 


B 2 (l + r 2 ) > 2r (^) > 0). Since r' —> 00 and p r —>- 0 (but still p <r < 0 and r' > 0), the first term in the bracket of 
Eq. ( [43] ) dominates unless B- 2 r = 0. 

Let us first assume that B 2 r = 0 at the pole r p . We then find 


B 


2 ,r 


= 0 

=> 

fi 2 = 

-/V P , 

(45) 

= 0 


(h = 

-#Rp, 

(46) 

= 0 


Ps = 

-/Vp, 

(47) 
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which leads to 


3H- 


= -f3 4 {r-r p )" 


(48) 


as well as 


B 2 = —/3 4 r p (r p - r ) 2 . (49) 

If r increases with time (which implies that dt > 0 since r' >0), then /3 4 has to be positive in order to get a positive 
H 2 . On the other hand, this would imply a negative B 2 and, thus, would violate the Higuchi bound. We therefore 
have to have a decrease of r with time (which implies contraction, dt < 0). Now, this is only compatible with negative 
values for /3 4 . Of course, such a model would be hard to believe in since it would contract at all times. But there is a 
more solid argument for ruling out these models: The contraction would lead to an increasing density. Since a root 
corresponds to a vanishing density, there must be a point of maximum density which always indicates a pole (see Eq. 
©)• Note, that we already excluded both r = 0 and r — ► oo as asymptotic states, which are the only ones that 
would be able to describe an infinitely large density. However, on this one, 7~L cannot vanish. Even though this second 
pole does not necessarily need to be an asymptotic point, H = 0 is required due to the bounce. Neither the positive 
nor the negative values for /3 4 lead to viable solutions and we conclude that 

7. Viability enforces a nonzero value for B 2iV around a pole. 

We are now allowed to assume H 2 r 7 ^ 0. Then the term —1 in Eq. (431 is negligible and, thus, 73 2>r has to be positive 
in order to get stability, i.e. w 2 < 0 . 

We will now study the expansion rate around the pole r p and check whether H 2 is positive. Note that (fh 2 ) r r 

does not automatically vanish since H 2 could become negative, too (which, however, would not correspond to physical 
solutions). However, the conditions for scalar stability (B 2 r I > 0), the existence of a pole (p r I =0with'H 2 | =0) 

and physicality (p\ r > 0 ) together with the assumption that (% 2 ) 
assume (H. 2 ) r \ r =0 which, together with H 2 \ r = 0 , implies 


7 ^ 0 lead to a contradiction. Therefore, let us 


1 


A = -3 r - 


-1 


( 2 /?! - /? 4 r p ) 


A = 3>V 


-2 


(Pi - 2/3 4 r p 


If we now assume that H . 2 is positive and nonzero at second order, then we need to have 

/ ti \ 2 

3 


(ylp rr2 ^ > (Pi + Pi rr l) = r P 3 (Pi + Pi r p ) (r - r p ) 2 + O ({r - r p ) 3 ) > 0 . 


However, this would imply that 


1 + r 2 
P,r = 2—^ 


(/3i + /3 4 r p ) (r - r p ) + O {{r - r p ) 2 ) 


(50) 

(51) 

(52) 

(53) 


only becomes negative when leaving the pole, if r decreases with time, i.e. in a contracting universe. We can now use 
the same argument that we used before and conclude that we need to reach a second pole which will either describe 
an asymptotic point or a bounce. Eq. (|52| provides the possibility of another point r P2 = —/3i/ {fi 4 r 2 ) at which the 
expansion stops but this cannot be a pole since then we would find 


P,r 


(/?i + /3 4 r p ) 2 (/?I + /3 4 r p ) 


7^0. 


/3i/3 ir% 

Our last chance are models in which "H 2 vanishes up to second order implying that 

Pi = ~Pi rl, 

P2 = Pat 2 


(54) 

(55) 

(56) 

(57) 


P3 = -pAl'p- 

These solutions lead to B 2 r = 0, which we already excluded earlier. Therefore, 

8 . A negative B 2 ,r around a pole leads to gradient instabilities whereas a positive value violates either the Higuchi 
bound or leads to unstable scalar perturbations. 


In combination with the requirement of a positive H 2 r . in order to get stable scalar perturbations this shows that 
every branch is plagued by either the Higuchi ghost or scalar gradient instabilities. 
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VI. CONCLUSIONS AND OUTLOOK 

We analyzed general models in singly coupled bimetric gravity around a FLRW background and found that all 
physical cosmological solutions that are not equivalent to ACDM have a period in time in which either linear scalar 
perturbations undergo a gradient instability or the Higuchi ghost appears. The condition for the absence of ghosts is 
surprisingly equivalent to r' > 0, which means that the ratio of the scale factors b and a has to increase as long as the 
Universe expands. Moreover, satisfying this bound ensures a positive lapse of f which is related to the absence of 
a helicity-2 ghost. 

In fact, all infinite branches suffer from the Higuchi ghost at all times and a ghost in the helicity-2 sector at early 
times, whereas in all finite branches, and even exotic branches that do not contain the limit r —> 0, there exists at 
least one epoch in which there is either a gradient instability in the scalar sector or a ghost appears. A schematic 
illustration of a typical phase space diagram with the forbidden regions is presented in Fig. [2j 

While the existence of a ghost renders the model unphysical and forces us to discard this type of model, unstable 
scalar modes will not necessarily rule out the theory. A Vainshtein screening may be able to prevent the scalar sector 
from getting unstable. Furthermore, this gradient instability is not present at all times. Every finite branch has a 
point in time at which the instability stops and the scalar perturbations begin to oscillate. As shown in [35], a small, 
but natural, Planck mass for f^ 6 would shift this gradient instability to very early times or even to energy scales 
above the cutoff of the effective field theory. In the latter case, the cosmological evolution would be very close to 
ACDM. On the other hand, if the instability ended between inflation and big bang nucleosynthesis, only very small 
scales would be affected |43J . These could, in principle, lead to a creation of many seeds for black holes. 7 

All models which we do not have to exclude due to the presence of a ghost will describe a phantom dark energy. 
Such a property would cause an anxious future in a ACDM model but not necessarily in bimetric theories due the 
time dependence of EOS corresponding to dark energy. In fact, it could cause welcome signatures that might allow 
observations to distinguish bimetric gravity from general relativity. 

Throughout this work we assumed a very simple, but well-motivated, type of bigravity. We considered a fluid 
that is only singly coupled to an observable metric and where both metrics are of FLRW type. Several extensions 
exist in the literature. One example would be the coupling of matter to both metrics and f^ v simultaneously 
EH S5MI, which, however, would introduce the BD ghost if the same matter sector is coupled to both metrics 
[1511501151] , Ghost-free (but not always with well-behaved cosmological solutions) scenarios exist if one assumes a 
coupling through a composite metric [10, 50] 150050] . But even the bimetric gravity with a standard matter coupling 
could allow for cosmological solutions without any gradient or ghost instabilities at the cost of giving up a FLRW 
background [Ml Ell. 
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